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Abstract. In the paper the rescaled occupation time fluctuation process of a certain empirical 
system is investigated. The system consists of particles evolving independently according to a 
symmetric Q-stable motion in R'', a < d < 2a. The particles split according to the binary critical 
branching law with intensity V > 0. We study how the limit behaviour of the fluctuations of the 
occupation time depends on the initial particle configuration. We obtain a functional central limit 
theorem for a vast class of infinitely divisible distributions. Our findings extend and put in a unified 
setting results of [S] and [22| . which previously seemed to be disconnected. The limit processes form 
a one dimensional family of long-range dependance centred Gaussian processes. 
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1. Introduction 

We consider a system of particles in R"' starting off at time t = from a certain distribution (to 
be described later). They evolve independently, moving according to a symmetric a-stable Levy 
process and undergoing the critical binary branching (i.e., or 2 particles with probability ^ in each 
case) at a rate V > 0. We concentrate on the case when a < d < 2a and the starting measure is 
an infinitely divisible point process M. We denote this system by N. It can be identified with the 
empirical process {Nt}^^^ which is measure- valued. For a Borel set A, the random variable Nt{A) 
indicates the number of particles of the system in set A at time t. We define the rescaled process of 
the fluctuations of the occupation time {XT(i)}(g[o i] 

(1) XTit):=jr£\Ns-X)ds, 

where Ft ■= T^^~~^^'^. The main object of our investigation is the limit of Xt as the time is 
accelerated, T +00. We are interested in how the limit depends on the starting measure M. 
In our paper we derive a simple integral criterion which determines the behaviour of the limit. 
The exposition will be clearer when put against the state-of-art in the field, consequently we start 
by recalling the most relevant results. The reader acquainted with them can proceed to the next 
paragraph. Studies of the fluctuations of the occupation time of branching particle systems started 



with pioneering papers |15| I16| . The field has received a lot of research attention recently due to 
the series of papers by Bojdecki at al. [71 El El [121 HOl El IH [H [H] and latter Milos [22l [26l 
[25l 123) . An analogous problem for a system of particles on the integer lattice was studied in [2]. 
The aforementioned papers concentrate on proving central limit theorems (usually in a functional 
setting) for various configurations of parameters (i.e. various values of d, a, branching laws and 
starting distributions). These findings were also complemented by corresponding large and moderate 
deviation principles e.g. [181 \T7\ [20l 124) . For a system with a critical finite variance branching law 
(e.g. the critical binary law) one can distinguish three main regimes of the limit behaviour: 

• low dimensions {d < a), when the system suffers from local extinction and direct study of 
the fiuctuations of the rescaled occupation time does not make sense, 

• intermediate dimensions {a < d < 2a), when the limit has a simple spatial structure 
(Lebesgue measure) and a complicated temporal one (with long-range dependence property), 

• large dimensions (2a < d), when the limit has a complicated spatial structure (5'(]R'^)-valued 
random field) and a simple temporal one (process with independent increments). 

The case of the intermediate dimensions, in which we concentrate in this paper, is interesting from 
a number of reasons. Let us recall results of [6l [22]. They considered two starting measures, a 
homogenous Poisson point process Poiss and the equilibrium measure Eq. The latter is the unique 
invariant distribution of the dynamics of the branching particle system. It exists under assumption 
d > a and can also be defined by (see |19) ) 

(2) Nt ^'^ Eq, t +00, 

where — ^"^ denotes the convergence in law and starts off from Poiss. Surprisingly, the limits 
of Xt are different in these cases (a so-called sub-fractional Brownian motion and a fractional 
Brownian motion, respectively). A natural question that arose was about the properties of starting 
distributions that are responsible for different forms of the limits. It became even more interesting in 
light of results of |14| (to be discussed in Remark 13. 4p . This paper brings an answer to this question 
for a vast and important class of starting distributions, viz. infinitely divisible point processes. 
We were able to factor out a single characteristic of a starting distribution M (named 'H{M) see 
(A2) in Section 12. Sp which determines the limit process. Precisely, its temporal structure (up to 
multiplicative constants) is described by a centred Gaussian process of the form 

where ^ is a sub-fractional Brownian motion, ^ is another centred Gaussian process, ^ and are 
independent. In this way we obtained a whole spectrum of limit processes. In our framework, 
Poiss and Eq starting distributions are special cases. Moreover, it is possible to give some intuitive 
description of T-L{M). An infinitely divisible point process is characterised by its canonical measure. 
Roughly speaking, this measure defines "clans of particles" which arrive to the system together 
in groups. The condition T-i{M) > implies that on average the size of a clan is infinite and 
asymptotically there are 7i{Q)r°^ particles in the ball of radius r. For the equilibrium measure these 
"clans" are exactly the clans which can be described in terms of genealogy bonds among particles 
(e.g. in |19| 130)). We do not know an exact explanation for such behaviour, although we expect 
that recurrence and persistence properties (e.g. like the one described in |29) ) play a crucial role. A 
clan has to be large enough so that its recurrence infiuences the limit. It is also worth mentioning 
that the limit process exhibits long-range dependence, which is also a consequence of the recurrence 



properties. Finally, we return to |14| : studying an analogous problem for particles systems without 
branching they obtained the same family of processes. We attempt to give an intuitive explanation 
of our results as well as understand the relation to the ones of |14) in Remark 13.21 and Remark 13.41 
The paper is organised as follows. In Section [2] we gather definitions and facts used in the subsequent 
sections. In Section [3] we present the main result of the paper, Theorem 13. H along with examples 
and comments. Sections U] and [5] are devoted to the proof of Theorem 13.11 Finally, in Section [6] we 
collected some basic facts about the equilibrium distribution. 

2. Preliminaries 

For the reader convenience we put a short index of symbols 
Symbols. 

A4{M.'^) - space of locally finite point measures on M*^ with topology of vague convergence 
Qx - Palm measure of random measure Q 
[fi] - set of atoms of a point measure fi 
A - Lebesgue measure 
9x - shift operator 

5(M'^), 5'(]R'^) - Schwartz space of rapidly decreasing functions and its dual space of distributions 

C([0, 1], <Y) - the space of ,^-valued continuous functions with the sup norm 

Pt{')-,Tt - density of the symmetric a-stable motion and its semigroup 

For a space of functions X by we denote the subspace of positive functions 

(•, •) - denotes either duality or integration 

~) ^ " use this notation when an equality or inequality holds with a constant c > which is 
irrelevant for calculations. E.g. f{x) ~ g(x) means that there exists a constant c > such that 
f{x) = cg{x). 

2.1. Functional setting. Process Xt is a measure- valued. However, from technical reasons, it 
is convenient to embed it in 5'(M'^) i.e. the space of tempered distributions dual to the space of 
rapidly decreasing functions 5(M'^). More precisely we will consider Xt to be a C([0, 1], 5'(M'^)) 
random variable and prove the convergence in law with respect to the topology of this space. 

2.2. Gaussian processes with long-range dependance. In our paper we will need two Gaussian 
processes. The first one is the sub-fractional Brownian motion. It was introduced in [5] where also 
its properties were studied. It is a centred Gaussian process with the covariance function 



(3) Ch{s, t) := s'' + t''-- [is + t)'' + \s-t 



h G (0,2). 



We will also need another centred Gaussian process with covariance 



(4) Chis,t) := sgn{h- 1) (s + t)'' - - t' 



h G (0,2). 



As indicated in |14) [after (2.4)] for h ^ 1 function is positive-definite since 




As the limit in Theorem 13. II we obtain a process which is sum of two independent processes as above, 
viz. it is a centred Gaussian process with the covariance 

e(s, t; h, A, B) = ACh{,s, t) + BChis, t), A>0,B>0, 

where h G (1,2). We already know that for A = l,i3 = 0itisa sub- fractional Brownian motion, 
moreover it is trivial to check that for A = 1,B = 1/2 it is a well-know fractional Brownian 
motion (see e.g. [28|). We denote by ^ a centred Gaussian process with the covariance function 
e. Such processes exhibit a long-range dependance which expresses the fact that the correlation of 
increments on distant intervals decays polynomially slow. We briefly summarise basic properties of 
^ in the following 

Fact 2.1. Let h G (1,2). Process ^ is self-similar, meaning 

{^iat),t > 0} {a^''^i{t),t > 0}, for each a > 0. 
Let us also denote the correlation of increments by 

R{u, V, s, t) := E{^u - ^v){Cs - Ct), <u <v < s <t, 

we have 

N i-r'"'^ ifB>0 
R{u,V,S + T,t + T)r. {^^^ ,^5^0, 

2.3. Symmetric a-stable motion. In this section we summarise analytic properties of a symmetric 
a-stable processes. This can be found in e.g.|3]. Let a G (0, 2] and {rjf}^ be the symmetric a-stable 
motion in starting from x. It is a homogenous Markov process with the transitions densities 
denoted by {pt}tyo- We also define the corresponding semigroup 

(5) Ttfix) := Ef{r,f) = [ pt{x - y)f{y)dy, t > 0, 

where / : h^R is a measurable function. We have also the following property 

(6) pt+s=Pt*Ps, 

where * denotes convolution (which is nothing else but the Markov property) . We will also need the 
following self-similarity property 

(7) pt{x) = t-^/y,{t-'/^x). 

The Fourier transform (i.e. f{z) = J^d e^^^ f{x)dx) of the semigroup Tis 

(8) 7I/(z) = e-*l^lV>). 

We recall that for any / G 5(M'^) we have / G 5(M'^). It is also well-known that for any / < a and 
any t >0 we have 

E\T]f\^ < +00. 



2.4. Infinite divisible point processes. The point processes theory is well-developed (e.g. |21|): 
in this section we excerpt facts needed in our work. We denote the space of locally finite point 
measures on R"' by and endow it with the topology of vague convergence. Our starting 
measure M is an infinitely divisible point process (IDPP). We assume that its expectation is finite 
and the intensity measure is A. For any measurable function tp : M!^ i— >■ M_|_ the following Levy- 
Khinchine-like formula holds 

(9) Eexp {- (M, V,)} = exp I - / (l - e"^'^''')) Q(d/x) I , 

where Q is a measure on A^(R'^) — {0} which is called the canonical measure and fulfils (Qvr^S 1 - e"^) < 
+00 for any Borel bounded set B and ttb '■ A4(M'^) — )• M is the mapping given by vrB(/i) = fJ-iB). 
We associate with Q a family of measures {Qx}x(^R'i called the Palm distribution which fulfils 

(10) / {fi,h)G{fi)Q{dfi)= [ h{x) [ G{fi)Qx{dfi)dx, 

for any /i : M'' — ^ M+, G : 7W(M'^) M+ which are measurable functions - see [211 12.3]. 
IDPPs may sound a bit abstract, but in fact they are not in the context of branching particle 
systems. We give now two examples playing a crucial role in our studies. The first is the Poisson 
point process (with intensity A), for which the distribution Qx is concentrated on 6x- The second is 
the equilibrium measure Eq defined by ([2]) , its Palm measure is described in Remark 13.11 where we 
will discuss these examples further. We will need shift operator 9x defined as 9xy := y — x, ?/ S M'^ 
for any B CM.'^ and measure /i G ^^(IR'^) we put 

exB:=B-x = {y-x:yeB}, 0,/i(-) := ^ o e-\-). 

We say that a point process is simple if the multiplicity of each atom is 1. A simple point process 
M is fully descried by the (random) set of its atoms [M], viz. 

[M] := |x € M'^ : 4 G m} . 

2.5. Assumptions. In this section we list assumptions required further and shortly discuss them. 
(Al) M is a simple translation invariant IDPP with the finite intensity measure A. We will denote 

it shortly by EAf = A or more precisely for any / E /^^(M*^) 

(11) E(M,/) = (A,/), E(Q,/) = (A,/). 

M has a finite second moment. We will assume a condition slightly stronger then usual, viz. 
for any f G 

E{M,ff < +00. 

For a measure M fulfilling (Al) we denote by Q its canonical measure. Q is also translation invariant. 
The family of the Palm distributions of Q will be denoted by Qx- The translation invariance of 
implies that Qx = dxQo- Moreover we denote the intensity measure of Qx by A^; := MQx 
(A2) Let {pt * Ao)(x) = f^jpt{x — y)Ao(dy). We assume that 

supt'^/"-^||pj *Ao|L < +00- 
t 



(A3) There exists a non- negative number denoted by 'H{M) or 7i{Q) such that 

i-^/"-! {pt, Ao) ^ niQ), as i ^ oo. 
(A4) There exist C, e > such that for any f G £^(]R'^) and x G R'^ we have 



Let us discuss shortly the above assumptions. We skip (A2) and (A4) which are technical and not 
very restrictive. (Al) is quite natural as A is preserved by the dynamics of the branching particle 
system (in expectation). In other papers in the field it is common to specify only a narrow class 
of the starting distributions (usually a Poisson random field or the equilibrium measure Eq). Our 
approach brings a novelty by admitting a larger class of distributions (the reader should be however 
warned that the infinite divisibility yields a poissonian structure lurking behind the scene - see \21\ 
6.2]). The detailed comparison to the known results is deferred to the next section. 
The assumption (A3) is the most important one. When l-i{Q) > it implies that particles appear 
at the beginning in "large clans". We postpone further discusstion of this phenomenon to Remark 
13.21 focusing now on giving some quantitive insights. The condition (A3) is roughly equivalent to 
the fact 

t~"Ao(-B(t)) converges, 

where B{t) is a ball of radius r. This equivalence is Tauberian in its spirit (though does not fall under 
the classical theory). For illustrative purposes we prove the following simple lemma, not pursuing 
analytical depth of the question 

Fact 2.2. Assume that there exist C > 0, e > such that 

(12) i-"Ao(S(t)) G {n{K) - Ct-^UiA) + Ct-'), Vt>i, 
then 

(13) t'^/"-^ (pt,Ao) ^ c„?^(A), ast^+oo, 
where = a pi{x)x'^~^dx. 



The proof is deferred to Appendix. 



3. Result 



Let us denote convergence in law by — t-"^ and recall the covariance functions from Section 12.21 We 
recall also the Gamma function T{z) = t^~^e~^dt. 

Below, we present a functional central limit theorem, which is the main result of the paper 

Theorem 3.1. Assume that a < d < 2a, and Ft = J'(3-<i/a)/2 _ N be a branching particle 
system starting off from a random point process M fulfilling assumptions (A1)-(A4)- Let Xt be its 
resettled occupation time fluctuations processes given by ([T]) . We have 

Xt^'^X, inC{[0,l],S'{R'^)) asT ^ +00. 

where X is a measure-valued centred Gaussian process with covariance function 

Cov{{X{s),^i) , {X{s),ip2)) = (A, (^i) (A, ^2) [KiChis, t) + n{M)K2Ch{s, t)] , 



where 'H(M) is defined as in (A3), /i := 3 — d/a G (1, 2) and 

r. ._ VT{2-h) _ 1 

1- 2"'-i7r'^/2ar(d/2)/i(/i - 1) ' ^' 

Remark 3.1. We now present a few examples. We keep the assumptions of Theorem 13.11 changing 
only starting distributions M 

• Let M be a Poisson random field. Assumptions (A1)-(A4) are easily verifiable and obviously 
T-L{M) = therefore limit process X has covariance function 

• Let M be a "compound Poisson random field with finite clans". It is a random point process 
obtained by replacing each atom of a Poisson point process by an independent copy of 
another point process. More precisely, let be a simple random point process such that 
EA^(R'^) = 1 and A^(]R'^) < C for a certain constant C > 0. Let P be a Poisson point process 
with intensity A. We construct M by 

a;G[P] 

where A^^' are independent copies of A^. It is straightforward to check that M is a simple, 
translation invariant IDPP with the canonical measure 

Q = f e^Ndx, 

where Af is the distribution of A^. As usual we denote the Palm distributions of Q by Qx- 
Our assumptions imply that intensity measures Ax := ^Qx fulfil Ax-(M'^) = 1. Using ([7]) we 
obtain 

Therefore assumptions (A2)-(A3) are fulfilled with HiQ) = 0. (A4) can be checked utilising 
condition A^(M'^) < C. Finally, using Theorem 13.11 we conclude that limit process X has 
covariance function 

Ki{X,ipi){X,ip2)Ch{s,t). 

which is the same as in the case of the Poisson point process. 

• Let M be the equilibrium measure defined by ([2]). It is known that it is a simple IDPP with 
a finite second moment; let us denote its canonical measure by Q. The Palm distributions 
Qx are equal to the law of 6x + , where is the limit of 

ef = Tc'^^^cds), 

Jo 

where i/ is a "Poisson clock" with intensity V, {Xg}^^^ a symmetric a-stable process starting 
at X (we think of it as moving backward in time), {Ct'^}t>o independent branching particle 
systems starting with one particle at y at time —s. is a clan of particles having a common 



(14) 



ancestor, for details and formalisation see [19^ I30|. We prove in Section [6] that the intensity 
measure of is [60 if irrelevant in limit) 

Ao(dy) = ^|yr%. 

Using ([7]) we can make explicit calculations (we denote Ci := 

t^/°-i(pi*Ao)(x) = Cit'^/"-^ / pt{y)\x-yr''dy = Cit~' [ pi{t~'/''y)\x - yr'^dy = {*). 

Now changing integration variable y t^^^^y we get 

(*) = Cit"'/"-^ / pi{y)\x - ti/"y|°-'^dy = Ci f pi(y)|t-^/°x - 
Using the last expression we easily verify (A2) and pass to the limit 

t'^'^-' (puAo) = * Ao)(0) ^ ^ / p^dy =: H{Q). 

^ jR<i \y\ 

We defer calculation of HiQ) to (j45p . Further, using Theorem 13.11 we obtain that the limit 
process has covariance function 

Ki (A,^i) (A, (^2) [Ch{s,t) + l/2ch{s,t)] . 

We have already noticed that Ch{s,t) + l/2ch{s,t) is the covariance function of the sub- 
fractional Brownian motion. 

A careful reader notices that we did not check (A4). Unfortunately, because of lack of 
manageable moments of order 1 + e it is very hard. Instead, using techniques of [22] one can 
check directly that converges to 0. This is the only place where (A4) is used. 

Remark 3.2. In this remark we would like to illustrate the theorem in an intuitive manner thus we 
do not pursue full mathematical rigour. 

(a) We were able to indicate parameter ■H(M) of the starting distribution, which is crucial for the 
behaviour of limiting process. Measures Q and describe particles which together enter the system 
at time t = 0. We will call such group of particles a clan. For rigorous treatment of clans in case of 
Eq see [30]. 

(b) Assumption (A3) gives a condition required to clans to have "influence to the limit". In order 
to explain it we refer to Fact 12.21 Under assumption 7i{A) equation (jl2p states that the (expected) 
number of particles in a clan has to be of the specified density. In particular this implies that the 
clan has to be infinite. 

(c) The limit process X could be written in the form X = ^A, where ^ is a real- valued centred 
Gaussian process with covariance function KiCh{s,t) + T-L{M)K2Ch{s,t). It has a simple spatial 
structure and a complicated temporal one. Let us recall that ^ is a long-range dependance process 
with the covariance decay rate described by Fact 12.11 The long range dependance stems from the 
recurrence properties of the movement. This could be considered at particles' level and at clans' 
level. The former can be explained as follows, a particle which occupies a given set at time ti is 
likely to "still be there (or come back)" at some later time t2- Of course the likelihood decreases 
with t2 — ti (our results indicate that the decay order is —d/a). The clans' level dependance is of 



slightly different mechanism. Again we consider a particle which visits a given set at time ti. It 
is accompanied by a large number of relatives from its clan (clans are infinite!). So even when it 
fails to be in the set at time t2 there is still a chance that some of its kinsman is there. This is 
reasoning is delicate, let us notice that the clans members evolve independently and the only source 
of dependance is through the fact that they arrived at t = in a cluster. Since the clan is infinite 
the family bonds are more strongly correlated than particles itself (i.e. when 'H{M) > the decay 
order is —d/a + 1). 

(d) The above explanation could be put in a slightly more formal setting (the reader might find it 
useful to read the proofs first and come back to this remark later). By the polarisation formula we 
have 

Cov{{XT{s),ip) , {XT{t),ip)) = [L{s, 0, T) + L{t, 0, T) - 2L{s, t, T)] /2, 

where L is given by (I36p . One is further referred to (i38l) . The term J3 capture the dependance at 
the particle level and results in Ch{s, t). The terms Ji, J2 describe the dependance on the clan level. 
To see we consider the "clan-less configuration", viz. the Poisson point process. One checks that 
Q = /jgd Sxdx and hence Ji + J2 = 0. J2 turns out to be negligible in limit. Finally we draw the 
reader attention to (i39l) which could be proved to converge to T-L{Q)ch{s,t). 

(e) Another way of explaining the additional term is to notice that {rij-, fJ-) in Ji is the same as 

r rTt2 i'Tt2 rTt2 

JR'' JTti JTti JTti 

where A'^'^ denotes the branching system starting from Nq = fi. Therefore {tlt, fJ-) could be considered 
as the occupation time of the mean density of the A^'^ (that is the mean density of a clan). In this 
interpretation, Ji + J2 is a variance, and consequently the term l-L{M)K2Ch{s,t) stems from the 
variability of the initial clan configurations. 

(f) We notice also that the limit process does not have stationary increments unless 'H(M) = 

We already know that this happens in case of Eq but now we see that it can happen more generally. 
We suspect that this may be related to the speed of convergence to Eq. 

Remark 3.3. The limit process X could be written in the form X = ^A, where ,^ is a real- valued 
centred Gaussian process with covariance function KiCh{s, t)+T-L{M)K2Ch{s, t). We observe a simple 
spatial structure and a complicated temporal one. Let us recall that ^ is a long-range dependance 
process with the covariance decay rate described by Fact 12.11 The form of the limit is intuitively 
explained with the notion of clans recurrence. In |29) properties of clans are described in the case 
of Eq measure. A clan here can be easily imagined once we consider an "infinitely old" branching 
particle system {A^tj^gj^. At time 0, the clans are maximal sets of particles having a common ancestor 
in past. For the system of Brownian particles in d = 3 (i.e. the only intermediate dimension for 
a = 2) any given ball is visited infinitely often by members of each clan. This results in smoothing 
of the spatial structure but introduces long-range dependance to the temporal one. The clans have 
not been studied in the case of a-stable particle systems and measures different then Eq but we 
believe that the picture is still valid. 



Remark 3.4. This research was inspired by developments of [14]. They consider the same problem 
for systems of particles moving in M and a class of starting measures defined by 

M = yy5, , 

jGZn=l ■> 

where 9j are independent copies of a non-negative integer-valued random variable 6 fulfilling < 
+00 and for each j € Z and k = 1,2,... we have a random vector p-j^ = (^pj, i, pI. 2^ ■ ■ ■ j Pk kj "^ith 
values in + l)'^ such that (9j, {p\)k=i,2,... ) are independent. They proved that for the system 

of particles starting from A4 and moving according to a symmetric a-stable motion without branching 
the limit process X is a centred Gaussian process with covariance function (Theorem 2.4) 

C (A, ^i) (A, ^2) HmChis, t) + 1/2 (Var0) Ch{s, t)] . 

This result is very similar to ours with the role of 'H(Q) played by Var^. From this point of view 
the interpretation given in Remark 13.21 (d) seems appealing. In |14| the system with branching is 
also studied but in this case the limit depends only on E^. As the authors explain in Remark 2.5 
"since the trees tend to extinction by criticality, the effect of the fluctuations of 6 become more and 
more negligible Our results could be seen as complementary to this fact and assumption (A3) 
is a condition which prevents the fluctuations to disappear. 

Remark 3.5. A natural question appears what happens when the limit in (A3) is infinite i.e. the 
clans are infinite and "denser" that (I12p . We conjecture that in such a case one needs norming larger 
than Ft = r(3-"'/")/2 and the hmit process X is a centred Gaussian process with covariance function 

C(A,99i) {X,ip2) Ch{s,t), 

where C > 0. That would indicate that the condition of assumption (A3) defines a point of a "phase 
transition" from the case of small clans and processes with covariance function Ch to the case of big 
clans and processes with covariance function Ch- 

Remark 3.6. [22] covers a wider class of branching law, namely critical branching laws with finite 
variance. We conjecture that our theorem holds in this setting (with possible change of constants) 
too. 

We also suspect that analogous results hold for systems with infinite variance branching law (like 
the ones in [81126]). 

4. Proofs 

We follow the path, which proved to be successful in solving similar problems. It stems from [6]. 
The description below is based on |25| Section 3.1, Section 3.2]. 

4.1. Scheme. To make the proof more clear we present a general scheme first and defer details 
to separated sections. Although the processes Xj- are signed-measure-valued it is convenient to 
regard them as processes with values in 5'(M'^). In this space one may employ a space-time method 
introduced by [3] which together with Mitoma's theorem constitute a powerful technique in proving 



weak, functional convergence. For a continuous, 5'(M'^)-valued process X = {Xt}^^Q we define an 
jS'(]R'^^^)-valued random variable X by 

(15) (^x,^y.= l\xt,'^{-,t))dt, 

we will refer to X as the space-time variable of X. 

Convergence. From now on we will denote by XtjX the space-time variables corresponding to 
Xt, X. To prove convergence of Xt we will use its Laplace functional 

Lri^) = Ee-(^'^'*), ^> e 5(M'^+i)+. 
For the limit process X denote 

L($) = Ee-(^'*) , $ G S{R'^+^)+. 
Once we establish the convergence 

(16) Lt($) ^ L{^), as T ^ +00, V$g5(Kd+i)^. 

we will obtain weak convergence Xt — X. Two technical remarks should be made here. We 
consider only non-negative The procedure how to extend the convergence to any <^> is explained 

in [6l Section 3.2]. Another issue is the fact that ^X^, is not non-negative (which is a usual 

condition to use the Laplace transform). The usage of the Laplace transform in this paper is 
justified by the Gaussian form of the limit. For more detailed explanation one can check [6l Section 
3.2]. Detailed calculations for this part of the scheme will be conducted in Section [5.11 
Tightness. By |3] we know that Xt -^"^ X together with tightness of {Xt^t ™plies convergence 
Xt ^'^ X in C([0, 1], 5'(M'^)). The tightness can be proved utilising the Mitoma theorem |27| . It 
states that tightness of {Xt}t with trajectories in C([0, 1], 5'(M'^)) is equivalent to tightness of 
{Xt,(p), in C([0,r],M) for every G 5(M'^). Using the simple trick (see [6][Lemma on p. 9]) we 
may also assume that (p > 0. Let us recall a classical criterion [1] Theorem 12.3], i.e. a process 
{{XT{t),ip)}^~^Q is tight if for any ti,t2 > and there exist constants C > 0, 7 > 1 

(17) mXrih), if) - {XT{t2), p)f < C{h - t2V. 
Calculations are deferred to Sections 15.21 

4.2. One-particle equation. In this section we will introduce an equation which is crucial for 
further analysis. The law of the occupation time for the system starting off from a single particle at 
X is described by 

(18) UM, (j;,r,t) := 1 -Eexpj-^ (A^f , ^' (•, r s)) dsj , * > 0, 

where denotes the empirical measure of the particle system with the initial condition A^q = 5x- 
More precisely A^^ is a system starting from one particle placed at x and dynamics described in the 
Introduction. The following lemma gives the announced equation 



Lemma 4.1. Assume that ^ G 5(M''+i)+ then 

(19) < Vvj, < 1, 

and V\i/ satisfies the equation 



(20) v^{x,r,t)= / Tt^s 







V n 

^' (•, r + t - s) (1 - t;* (•, r + t - s, s)) - —v^ {■,r + t-s, s) 



(x) ds. 



The proof follows by use of the Feynman-Kac formula and can be find in [6l Section 3.3]. is 
quite cumbersome to deal with hence we approximate it with n,j, defined by 

(21) n^{x,r,t) := [ 7^_,^'(-, r + t - s)ds, ^ G cS(M"'+^), x G M"', r, t > 0. 

Jo 

Intuitively, nijr was obtained by dropping quadratic terms in (I20p as they do not play a role when ^ 
is small. 

Notation. From now on we fix non- negative $ G 5(M'^'''^)+ and prove convergences announced in the 
scheme in Section 14.11 To make the proof shorter we will consider $ of a special form. The proof 
for general ^ goes the same lines but is more notationally cumbersome 

(22) $(x, s) := ^pix)ij{s), if G 5(M'^)+, V G 5(M)+. 
We also denote 

(23) (fT (x) := -^(p (x) , x{s) ■■= 'ilj{u)du, XT ■= X (^j; 
We write 

^'(j;,s) := ip{x)x{s), 

(24) *T {x, s) := {x, ^) = ^t{x)xt{s). 

Note that ^ and are positive functions. In the sequel, we also write 

(25) UT(a;, r, t) := Uijiy (x, r, t) and vt{x) := vt{x,Q,T), 

and we define in the same way nT{x,r,t),n'r{x). Finally, by (j20p . (f2T]) for ^' > and by the fact 
that the potential operator of a-stable motion exists for d > a, it is straightforward to check 

(26) < vrix) < nrix) < V^gjjd, 
and some C\ii > 0. We will also use the following simple estimate 

1 

(27) nrix) - vt{x) < —nxix) + TT-snT{x,T - s, s)ds. 

Jo 

which again follows from ([20]), dH]) and (l26|). Let us fix ^ G 5(1^^^)+ and denote v{6) = vgq,. 
Obviously v{0) = 0. In the sequel we will need derivatives of v with respect to 6. Using (f20|) it is 
easy to calculate that (we omit arguments and integration variables) 

(28) v'{0)= I Tt-s[^{i-v{e))-e-^v'{e)-Vv{e)v'{e)\ds. 



When ^ = we have 

rt 



(29) v'{0){x,r,t) = [ Tt-s'^{x,r + t-s)ds = nq,{x,r,t). 

Jo 

We win also need the second derivative 

v"{d) = [ Tt-s [-2^v'{e) - eW{9) - Vv{e)v"{d) - v{v'{e)f] ds. 

Jo 

When ^ = we have 

(30) v"{0){x, r,t) = - [ Tt-s [2^'(-, r + t- s)n^{-, r+t-s,s)+ Vnq,{-, r + t-s, sf] ds. 

Jo 

By psp and properties of the Laplace transform we have 

E {N!, ^ (., s)) ds^ = -v"mx, 0, t). 

Formally the equation is true for ^ € 5(M'^'^^) + . Exploiting the Lebesgue dominated theorem one 
can easily argue that it is also true for ^ of the form ^{x, s) = (p{x)l[i^^^^^ yielding 



t2 



2 



(31) E(^y^ {N^,ip)dsj = W(0)(x,0,t), 

with v" defined according to (|30p with this special ^. 

4.3. Miscellaneous. In this section we gather some simple facts used in the proof below. All 
functions in this section as assumed to be measurable 

Lemma 4.1. Let us assume (A2). There exists a constant C > such that for any ip : M.'^ i-^ M_|_ 
we have 



supt-^/^-i / {Tt^,fi)Qoidfi) < cM^. 

Proof. Using ^ we have 

/ {Tnp,n)Qoidfi) = Tnp{x) Ao{dx) = / pt{x - y)ip{y)dyAo{dx) 
Jm(R'') Jr'' Jr'' Jr'' 



ipiy)ipt * Ao)(dy) < \\pt * AoHoollv^lli- 

Now the proof follows directly by (A2). □ 
Lemma 4.2. Let us assume (A2) and ut he given by (j25p then 

{nr^n) ^0 in C^{Qx)- 
Moreover the convergence is uniform in x. 



Proof. Using similar calculations as in the previous lemma we have 

{nT,fJ,)= / {nT,fJ.)QxifJ') = / nT{y)A^{dy) = / / TsipTiy)<isAx{dy) 

[ Ts^Tix + y)dsAo{dy) = F^^ [ [ <p{x + y){ps * Ao){dy)ds < [ s^-^ds 
Jo Jo Jr-^ Jo 

□ 

Lemma 4.3. Let fi, /2 : M'^ M+ . Then 

/ / Ts^h{x) {Ts^f 2, fj)Q x{d^x)dx = I (Tli+^aF*,^) Qo(d/i), 

ly/iere := Jj^d fi{-z + 2;2)/2(2:2)d2;2. 

Proof. We recall that Q^,' = GxQo^ hence 

/ / Tsjiix) {Ts2f2,fJ')Qxid.n)dx = / Tsjiix) {Ts2f2,OxfJ')Qoidn)dx = {*). 
Jr^ JMiR'') JR"^ Jm(R'') 

Using ([5]), the fact that Qo is a simple point random measure and the Fubini theorem again we 
obtain 

(*) = / PsT_{x - zi)fi{zi)ps2{x + y - Z2)f2iz2)dxdzidz2Qo{<ifJ') = {*)■ 

Jm{R'') Jr'' Jr-^ Jr'^ 

In the next step we substitute x ^ x + zi, use the semigroup property ([6]) and symmetry of Ps to 
obtain 

(*) = / Ps^{x)fi{zi)ps2{x + y + zi-Z2)f2iz2)<ixdzidz2Qo{dn) 

JMiR-^) ^gj^] JR'' JR'^ JR"^ 



Further, substitute zi zi + Z2 and we write 



/ / fi{zi)psj^+s2{x + y + zi - Z2)f2{z2)dxdzidz2Qo{dlJ,) 
TTi JR'' JR'' 



(*) = / y2 fiizi + Z2)p(si+S2)iy + zi)f2{z2)dzidz2Qo{dn) 

JAI(Md) ^^^^^^ jRd 

= / ^2 Pisi+S2)iy + zi)F*{-zi)dziQo{dfi) = / (T^i+sai^*,^) '3o(dM)- 
Jm{R'^) ^gj^j JR'* JA1(IR'*) 

□ 

We present also a simple fact concerning moments of a IDPP random measures. 



Lemma 4.4. Let M be a IDPP random measure with canonical measure Q. Let f -.W he 
such that E (M, /)^ < +00 then 



E{M,f)= / (/,/u)Q(m). 



2 / 

\2 



\Jm{M'') J \Jm(R'') ) 

We omit the proof which fohows by standard appHcation of the properties of the Laplace transform. 

5. Calculations 

In this section we execute calculations announced in the proof scheme described in Section 14.11 

5.1. The Laplace transform. In this section we calculate the Laplace transform of the space-time 
variable corresponding to X^- Let us recall definitions ([T|) and (jlSp . Performing simple calculations 
we obtain 

L($) = Eexp {- (It, ^)} = exp (^t(-, s), A) dsj Eexp |- (iV„ ^'t(-, s)) dsj. 



Next we condition with respect to the starting measure M and apply ([18 



Efexpj- / (iV„x&T(-,s))d4 hl\= n Eexpj- / (iV,^ I/tO, s)) ds) = TT (l-^r(^)) 
Now using identity OxgIM] ~ vt{x)) = exp {(M, ln(l — vt))} and Q we obtain 

L($) = exp I / ^t(x, s)dsdxl exp | - / (l - e^Hi~VT),,.)\ g^^^) I _ 

UR'* Jo J [ Jm{M'^) ^ ^ J 

Let us denote the exponent of the expression above by /(T), viz. Eexp |— (^Xt, *^^| — ^xp {/(T)}. 
We introduce the following decomposition 

I{T) = / / "^t{x, s)dsdx — / ?;T(2;)d2;+ / ln{l — vt{x)) + VT{x)dx 



Ii{T) l2{T) 

I 



[ (1 + (ln(l - VT),fi) - e<i°(i-''-)'^>) Q(dM), 



where between the lines we used the identity {fi = j^d f{x)dx (A is the intensity measure 

of Q). The proof will be completed once we obtain 

Ii{T)^—[ [ [ [ ^{x,t)^{x,t)dxdyCh{s,t)dsdt, as T — 7- 00. 
2 Jo Jo Jr'^ Jm.'^ 



hiT) ^0, as T ^ oo. 

(32) h{T) -HiQ)— f [ [ [ ^{x,t)^{x,t)dxdych{s,t)dsdt, as T ^ oo. 
where Ki and K2 are the same as in Theorem 13.11 For now using the above hmits we deduce 

I{T)^\C C f [ ^x,t)^x,t)dxdy[KiCh{s,t) + n{Q)K2Ch{s,t)]dsdt. 
^ Jo Jo JR'i JR'i 



It is easy to check that this concludes the proof of (I16p . Therefore according to scheme from Section 
14.11 we are left of the proof of tightness. This is deferred to Section 15.21 and now we proceed to the 
proof of the converges announced above. 

Convergence of Ii. This convergence can be recovered easily by inspecting [6l (3.23)-(3.30)]. 
Convergence of 12- Inequality (f26|) entitles us to use of the inequality | ln(l — a;)+a:;| < x'^,x G (0, 1/2). 
Whence by (f26|) and ([2T|) we have 



1^2 (T) I < / v^ix)dx< / n^(x)dx 



j ^ (^j Ts(pT{x)ds^ dx. 



Using (l8|), rearranging, using the Fubini theorem and recalling (p3|) we arrive at 

|/2(T)|< / r r e~^^\^\"e-^-\^\"ds,ds2dz<F-' f Mg! (1 _ g-^l^r)' dz. 

Jo Jo Jk" kr" ^ ^ 

We now introduce an inequality which will be used a few times. Namely for any 6 € (0, 1] we have 
(33) 1 - e-" < x\ V,>o. 

Now using the inequality with 5 = 1/2 and recalling that = T^~'^/^ we obtain 



1/2 (T) I < r'^/"-2 I 

Jw 



-dz as r — )■ +00. 



where the last integral is finite since d > a and ^ € 5( 
Convergence of I^. Using the Palm formula (llOp we write 

r m /■ Wi ( n f l + (ln(l-z;T),/i)-e<''^(^-"^)'^) ^ , , , . 
h{T)= ln(l -t;T(x,0,r)) / — — ^ — ^ Q^(d/i)dx, 

where Qx is the Palm measure of measure Q. Let us recall that nT ~ vt^ we know also that for 
small X we have ln(l — x) ~ —x and ~ ~f ■ Hence I3 is well-approximated by 

hi{T) := -\ I nrix) I {nT{x), fi) Qx{dfj,)dx. 



2 



First, we will prove first the convergence of /31 to the same limit as in (|32p . In the second step we 
win prove that |/3(r) - /3i(r)| 0. We apply definition of ut (recall ([21]) and ([25])) 

hi{T) = -\i / Ts^LpT{x)xT{si)dsi i (/ Ts2'^TXT{s2)ds2,li) Qx{dii)dx. 



Using the Fubini theorem we get (ah functions are positive) 



hi[T) 



Xt{si)xt{s2) 



'0 JO 

Let us substitute si — )• Tsi S2 — )• Ts2 and recah (|23p . Using them we acquire 

hi{T) = -\'^[ f x(.si)x{s2) [ [ TtsM^) {'TTs2'P,f^)Qx{dn)dxds2ds 

^ Jo Jo JR'' JMiW^) 

Using Lemma [4.31 with fi=f2 = ^ and := + Z2)(p{z2)dz2 we obtain 

hi{T) = -]-^( [ x(.si)x{s2) [ {Ts,+s2F*,n)Qo{dfi)ds2dsi. 
^ j't Jo Jo jm(m.'^) 



T 

We also recah that F| = T^"'^/", hence 



hiiT) 



1 rl 



Jo 



X(si)x(52)(si + S2)'-''/"(r(si+S2))''/"-' 



(T, 



By (A3) the inner expression converges to "^((5)11^*11^ = \\(p\\i and (A2) vindicates use of the 
Lebesgue dominated converge theorem, therefore 

hm I.,,{T) = -l-n{Q)M\i C l\{si)x{si){si + S2f-'"^ds2dsi. 
Let us now consider the inner integral and recall ( I23p 



1 rl /•! /•! 



/ / x(5i)x(si)(5i + S2)'-''/"ds2dsi= [ [ [ [ V(ni)V'(n2)dnidn2(si + S2)'-''/°ds2dsi 

Jo Jo Jo JO J Si J S2 

l[si,l](^^l)l[s2,l]('"2)^('"l)V'(^^2)(si + S2)^~''^°duidti2ds2dsi 
"1 rl 

l[si,l]('"l)l[s2,l](^2)(si + S2)^"'^/"ds2dsidliid'U2 



JO JO JO 
1 rl 

lp{ui)lp{u2) 



Jo 



7' 

Jo 



1 /•! 



JO 



(2-d/a)-W / ^K)^(^Z2) / l[si,i]{ui) (si + n2)2-'^/ 



\2~d/a .P'-d/a 



dsiduidti2 



{{'i-d/a){2-d/a)y^ 



1 rl 



Jo 



V'(ni)V'(n2) 



[Ui + U2 



duidti2 



1 rl 



JO 



{{3 - d/a){2 - d/a)) M / ■ip{ui)ip{u2)ch{ui,U2)duidu2, 



where h = 3 — d/a. This is exactly the same as limit (I32p . To conclude we still have to prove that 
|/3(r) — /3i(T)| 0. This will be done in a few steps in which we will gradually change formula 
from /3 toward /31 proving each time that they are asymptotically the same. Let us denote 

h2{T):= I Hl-VT{x))f -]-{\n{l-VT{x)),^l)Qx{d^J)dx. 

jRd. Jm(R'^) ^ 



Therefore we may estimate 

(34) h{T)-h2{T)= f H1-vt{x,0,T)) [ ZTQ,W)dx, 



where 

(ln(l - 7;t),/") 2 
Let us notice that (|26p and Taylor expansion imply that \Zt\ ^ min((?i5", /i) , (n-r,^)^), therefore 
by Lemma 14.21 we know that {nT,fJ') — >• in probability under Qx uniformly in x. By the diagonal 
procedure we may construct Ht — ?• such that sup^. QxiinT, fJ-) > Ht) — )• 0. We now write 

|/32(T) - /3(r)| = / ln(l - Vt{x)) [ ZTl{(„^,^)<H^}Q..(d^)dx+ 

/ ln(l - 7;t(2;)) / ZtI 

jRd Jm(R'') 



We have | ln(l — x)| < Using this the above considerations, (j26p and the Holder inequality we 
get 

|l32(r) - l3(r)| <Ht [ ln(l - vt{x)) [ {nT,fi) Qx{dfi)dx 

+ suvQx{Zt > HtY'^^^''^ [ ln(l - vt{x)) \ I Z^+'Qx{dfi) 
X jRd yJM{R<^) 

where e is the same as in assumption (A4). Applying (A4) and | ln(l — vt)\ < we arrive at 
|/32(r) - h{T)\ < Ft/31 (T) + supQ,.(Zt > HTY''^^+'h^i{T) ^ 0. 

X 

Next we define 

hz{T) -.= 1 [ VTix) [ (ln(l - vt), iA Qx{d^^)dx. 



1/(1+^) 

dx. 



2 

Using the inequality | ln(l — x) + x\ ^ x'^ valid for x £ (0, 1/2) in conjunction with (|26p we get 



|/33(T) - /32(T)| < / VT{xf f (ln(l - VT),fi) Qx{dfi)dx. 

JR'' Jm{R'') 

Applying (|26p and calculating as before we have 



|/33(T)-/32(r)|<^^^0. 

hi{T) := -]- I vt{x) I {vT,n) Qx{dfi)d^ 
^ JRd JA4(R<*) 



We define also 



By (p6|) and the Taylor formula we have | log(l — vt{x)) + VTix)\ < vt{x)'^ < Frp^nT{x). Now 
analogously as before we prove 

|l33(r)-/34(T)| ^0. 



Further we put 



2 



Using inequalities ([27|) and ([26]) it is easy to get |/34(T) — 135 (T)| < /35i(r) + 1352 (T) where 

/35i(r) = / ^nT(x) / {vT,fl) Qx(d/i)dx < ^/3i(r) ^ 0. 
and in the second term is (we use ([26|) under integral /_;vi(]gd)) 

/352(T) < / r Tt^shUx, T-S, S)ds I 



(nT,/i) (5x-(d/i)d3;. 



It is easy to check that nj'{x,T — s,s) < nxix). Using this, (j2ip and the Fubini theorem we write 

h52{T) < / / / / 7T-Sl^^T(2;) (T^aV'T,^) Qx(d/u)dxdsids2. 
JO JO JR'i Jm{R'') 

Using Lemma [431 with /i = n^, /2 = ifr and F^{z) = J^d n^i—z + Z2)¥'T(-Z2)dz2 we obtain 

(35) h^2{T)<r ri {TT-s,+s,F^,^i)Q^{di^)dsids2. 

Jo Jo Jm{K'') 

By Lemma [4. II we have 

/352(T)< r r{T-s + u)'-''/''\\F^,dsds2Qo{dfi). 
Jo Jo 

Using simple calculations and (j8)) we get 

II-^tIIi = / nT{x)'^dx / ipT{x)dx ^ F^"^ / / / %j.i^{^)%i2'Pi^)^''^i'^''^2'ix 



2 



1 

Let (|33p with <^=l~|^ + £- One easily checks that for e small enough we have 6 S (0, 1). Hence 
by ()33p follows 



|i?*|| < _p-32i(2-d/a)+2e /" ^ j^(d/a-5)/2+2e 

I Till T I \ ^\d—2ae 

Jw^ \Z\ 



Further 

rT rT 



h52{T) < rW"-5)/2+2. f f (^T-Si + 52)'-'/"dsids2Qo(dM) 

JO JO 

^ rp{d/a-5)/2-2ej^3-d/a _ j^{l~d/a)/2+2t _^ q 

if only e is small enough. Therefore 

|/34(r)-/35(T)| ^0. 



Finally, let us notice that by inequality (|27|) we have |/35(T) — /3i(T)| < l36i{T) + /362(r) where 
Imi{T)= f nrix) [ ( ^nrix), fi) Q,{dfi)dx = 2^^^ ^ 0, 



and 

rT 



h&2{T) = I nrix) / ( / TT-snT{-,T - s, s)ds, fi) Qx{dfi)dx. 
Using the Fubini theorem, inequality n'^{x,T — s,s) < n'^{x) and (|2ip we get 

h62{T) < / / / / TsiVT{x){TT-s2nT,fJ')Qx{dfi)dxdsids2. 
Jo Jo Jr'' Jm{R'') 

Using Lemma with /i = ifx, H = and F^{z) = J^d n^i—z + Z2)^PTiz2)dz2 we obtain 

l362(r) < r r [ {TT^s,+s,F^,f^)Qo{dfl)dsids2. 
Jo Jo J7M(K'*) 

This is the same as (|35p . Therefore /362(r) ^ and the proof is concluded. 

5.2. Tightness. In the whole section we fix < ti < i2 ^ 1- Our aim is to compute the second 
moment of the increment, viz. E((Xx(t2)) 9^) — i^Ti'ti),^))'^ and proof pT|) . In this section we use 
the following definition of ht 

rTt2 

nxix) = / Tu(pT{x)du. 
Jrti 

Obviously we have 

(36) L{ti,t2,T) := E {{XT{t2),^) - {Xriti), = E (^j^'' (iV„ - A, ^t) di^j 

= E (^^^ ' (iV„, c^t) du - {Tt2 - Tti) (A, (^T>) 
= e( f \Nu,y^T)du) -2{Tt2-Tti){X,^T)^( f \N.,,ipT)) du + {{Tt2 - Th) {X,y,T)? . 



\JTti / \JTti J 

It is easy to prove that A is the intensity measure of A^^^ (which follows by the fact that A is the 
invariant measure of the a-stable motion). Using the Fubini theorem and equation (|lip we obtain 

E ^ /j^*^ (-^«; V't) du^ = J^^*^ E (A'u, v^t) du = {Tt2 — Tti) (A, ifx)- Therefore we have 
(37) L{ti,t2,T):=E(^£\Nu,VT)du^ - {{Tt - Tt) {\,^T)f . 



Let us deal with the first term; we first condition with respect to the starting measure M 
Li{h,t2,T) := E (^j^^ ' {Nu, Vt) dv^j = EE ^ (^J^^ ' {N^, ^t) dn^ 



2 




) 





EE 



where by we denote a system starting from 5x- Let us calculate (slightly abusing notation we 
use E also to E(-|M)) 



Li{tut2,T)=E V E / 



xe[M] 



Tt2 S 2 

(iV^,V?r)dn 



+ E ^ E 



Tt2 



Tti 



Tt2 



Tti 



{Ny,^T)du]. 



The systems and are independent as long as j; 7^ y, therefore 



E 



Tt2 



{K,ipT) du 



Tt2 



Tt2 



{Nl ^T)du\=E{ / ^T)du\E{ / (iV,^ ^t) du = (*). 



Tt2 



TU 



We know that E (A'^^, (fx) = Tu'-PT{y) (this is standard and can be found e.g. in |19]). Due to it 



(*) 



Tt2 



Tu^T{x)du 



Tti 



Tt2 



Tti 



Tu^T{y)du = nT{x)nT{y) 



Now the second summand of Li can be recalculated as (by assumption (Al) this expression is finite) 
E V nT(x)nr(y) - V {nT{x)f = / (/x, nr(a;))^ - (/x, nr(a;)) Af (d/i). 



/EG[M],ye[M] x&[M] 

By Lemma [4.41 and equation (jlip this is equal to 



/ (;U,nT(x))^(5(d/i) + ( / nT(2;)dx I - / n\{x)dx. 
We also use ([HI]) with ^'(x,s) = (y9'p(j;)l[j.j^ j^f^]!*) ^^d (fTT|) again in order to get 
Li{ti,t2,T):= [ -v"{0){x,0,T)dx+ [ {^x,nT{x)f Q{d^l) + ( f nT{x)dx] - [ 
Let us notice that J^d nT{x)dx = {Tti — Tt2) (A, (fr), hence L given by ([37|) writes as 
(38) L{ti,t2,T) := I {nT,fifQ{dn)- I n'^{x)dx + [ -v"{0){x,0,T)dx . 



ri^{x)dx. 



MT) 



MT) 



MT) 



We start with a simple case of J2. We have 

fTt2 \ 2 pTt2 rTt2 



/> / />_Z t2 \' rit2 /•lt2 r 

J2{T) = I ( / Ts'^T{x)ds\ dx = / / TsjipT{x)Ts2^T{x)dxdsids2. 

JR'' XJTti J JTti JTti JR'' 



Changing variables si — > Tsi and S2 Ts2 and using (j23p we get 

'■t2 (-12 



rt2 rt2 r 

J2(r)=T^Fy^/ / / TTsi^{x)Ts2^ix)dxdsids2 
Jti Jti JR'i 



Applying the Fourier transform ([8]) we obtain 



Jti Jti JR'i 



Further 



We use inequality (I33p with 5 = (3 — (i/a)/2 



J2(T) < (i - / l^dz < (i - 



Now we deal with Ji. Using the Palm formula (|10p we obtain 

Ji(r) = / / {nT,fi)Qx{'ifJ.). 

JR'' J Mi 



Using the definition of tit we write 

rTt2 r I rTt2 



Jl{T) = / / Tsi'-PT{x)dsi / ( / Ts2^T{x)ds2,M Qxi'^lA'^x. 

JR'iJTti J MiR"^) \JTti I 

Now we apply the Fubini theorem and recall (f23|) to get 

(■Tt2 rTt2 r r 

Ji{T) = F^^ / / TsMx) {Ts2^{x),li)QMlj)'^xdsids2. 

JTti JTti JR'' Jm{R'^) 

Using Lemma [O] with fi = ip^ f2 = ip and F*{zi) = J^a fi-^i + Z2)^iz2)dz2 we get 

/•Tt2 rTt2 r 

(39) J,{T) = F^^ / / {rs,+s2F*{x),fi)Qo{dfi)dsids2. 

JTti JTti JMiR'i) 



Now by Lemma [4.11 we have (we recall that h = {3 — d/a)) 



i-Tt2 rTt2 
JTti JTti 

= (^(jt^ + Tt2f - {Tti + Tt2f - {Tt2 + Ttif + {Tti + Ttif'^ 

= {2t2f - 2(il + t2f + {2hf < {t2 - hf. 



Finally, we deal with J3. Using equation (j30p we write 

J^{T) = 2f [ TT^s[^{-,T-s)n^,{-,T-s,s) + Vn^{-,T-s,sf]dsdx 

JRd Jo 

^ / "^{x,T — s)n^{x,T — s, s)dsdx + / / n^{x,T — s, s)'^dsdx . 
JR'i Jo Jw^ Jo 

J3l{T) J32(T) 

We start with J31, let us recall that = (/3T(2;)l[Tti,Tt2] ('^))' substitute s ^ sT and ^ Tit 

and apply (|23]) 

hi{T)=l / / v5T(a;)l[Tti,Tt2](s)'7^-n'/?T(2;)l[Tti,Tt2](^ - ^)d^^dsdx = 

JM'' JO Jo 

[ [ f <^(a^)l[ti,t2](^)'7T(s-«)y'(a^)l[ti,t2](l -'")d^^dsdx. 
Applying the Fourier transform ([8]) we get 

JsiiT) = ^ [ f r |^(z)pl[,,,,](.)e-^(-")W"l[,,,,,](l - ^)d^d.dz. 

i^T JK'* JO JO 
We may easily estimate 

^31 (T) ^l[..,.2](-)e-^(^-*^^l^l" (1 - e-^(*-*^)l^li d.d. 

i'T JK"* JO 1^1 ^ ^ 



T 

we use inequality (j33p in order to get 

T 



l+<5 



as one can easily find a proper 6 > 0. Let us now proceed to J32. Using (I2ip and changing variables 
s Ts, u Tu 



h2{T) = j^J^ Ts-u^T{x)l[Tt,,Tt2]iT - U)du^ dsdx 




1 / rs \2 

TT(s~uMx)l[t,,t2]i'^ -u)du] dsdx. 



Using equation (jH]) we obtain 

-^32^^!^ / /'|^(z)|2(re-^(^-")W"l[,^,,](l-n)d..)'d.dx. 

Next, we write 

T^3 rl rs ps p 
< — / / / p-T{2s-m-U2)\z\- . 



«^32(T)<-2/ / / / e-^^^^-^^-"^^l^l l[i,,i,](l-ni)l[,,,i,](l-n2)dnid^.2d.dz. 

Jo Jo Jo JR'i 



T 

Integrating with respect to z yields 

"1 fS rs 



•^32(^)^ /„ /„ j {{2s-Ui-U2)r^'''\t,M^l-Ui)l^t^^t^-^{l-U2)duidu2ds, 

Further 



^0 JO 



rl-ti rui rl 

J32{T)< / / / {2S - Ui - U2)~''/''dsduidu2 ^ 

/ {2-ui- U2)-'^'''^^du2dui - / {ui- U2Y'^''^^^du2dui. 



Let us denote S = t2 — ti and write 

7"!— i2+t^ pui 
1 

It is easy to check that by rHopital's rule we have 
/(<^) ~ lim i I 

5^0+ 



1 t2+o rui 

m ■■= 7^ I / (2 - «i - U2)-^/°+Mnidn2. 

t2 Jl-t2 



1 / rl-t2+5 \ 
lim f{5) ~ lim — — / {i + t2-5- U2)-^^"^'du2 < c. 



Hence we get Jf^~^^ JTi^i^ — ui — U2) "'^"'''^dnidM2 < (t2 — ^i)^ ■ The other term is estimated by 
(^xi - n2)i-'^/°dn2dni ^ {m - tif~'''''dui ~ (t2 - ti)^'"''/". 

-J. - i J ti 



6. Equilibrium measure and its Palm measure 

Let us recall the description of clan decomposition of the equilibrium measure given in Remark 
13.11 In this section we will compute the Laplace transform and the first moment of , see equation 
(fT4l) . To this end we first consider the following Laplace transform 

w{x,t) = l-Eexp{-{N^J)}, 



(40) 



where A^^ is the branching particle system with starting condition Nq = 5x and / : M i-^- is 
some measurable function. It is known |19| (2.3)] that 

(41) w{x, t) = Tt{l - e-^(-))(x) - ^ r Tt-sw{x, sfds. 

^ Jo 

Using this equation we will now calculate the Laplace transform of ^f. Let us denote by Poiss a 
realisation of v 



(42) 



H{t, x) := Eexp {- (Cf, /)} = EE[(exp {- (^f , f)})\ Poiss, X' 



Further, let ti,t2, ■ ■ ■ ,tn be the times when Poiss ticked (n is an a.s. finite random variable ). We 
have 



H{t,x) = EE 



exp 



\i=l 



Poiss, X"" 



E 



nE[(exp{-(c;:"'-,/)} 



Poiss, X' 



1=1 



Utilising (|4U|) and the definition ([S]) we arrive at 

H{t,x)=E(jlE{{l-w{Xt^,ti))\ Poiss, X"")^ =E(j[ru{l-w{x,U))j . 

Using the properties of the Poisson random fields we obtain 

(43) H{t,x) =expl^-V J^Tsw{x,s)ds^ . 

Expectation of We denote w{x,t, 6) given by (j40p defined for 9f. Obviously, it fulfils an analogue 
of dH]). We easily check that 

w'{x, t, 6) = 7I(/(-)e-^^(-))(x) -V f' Tt^s [w'i; s, e)w{; s, 9)] {x)ds, 

Jo 

w'ix,t,0) = Ttf{x). 

We can calculate the moments of (^f,/) by differentiating H{x,t,6) which is an analogue of ()^2|) 
defined for Of. In fact we may write 

H{t,x,9) :=exp|-yy" Tsw{x, s,e)ds^ . 

Differentiation with respect to 6 yields 

H'{t,x,9)= (^-y Tsw'{x,s,9)d^ exp|-y^ Tsw{x, s,e)ds^ . 



which evaluated at = is 

H'{t,x,0) = -V I T2sf{x)ds. 
Jo 

Using the properties of the Lebesgue transform we read 

mfJ) = V fT2sf{x)ds 
Jo 

Applying the monotone convergence theorem and using the fact that system is persistent (see [19]) 
obtain 

(44) E(e^,/> = ^ Pr./(x)d. = ^ / ^^dy, 

where c„,rf = T (^) (2°7r^/2r(a/2))-i (see ^ (1.32)]). 

Now, we would like to calculate f^^ Pi{y) \y\d-a dy. Let us denote 

= p,{y)—J——dy. 

jRd. \x — y\ 

F(^) = c-i.-N"J_. 



Using dS]) and (l44|l we easily obtain 



Using the inverse transform we conclude 



"^^'^ (2vr)'^yK/^^ i2nrL kl" (2vr)^ r(d/2) io 



Now we substitute r — ?> s^/" getting 

'?r~^ ^dl2 /■+oo i 

(45) F(0) = ^. / e-^l.(^-)/".V«-id. 

(27r)'^ r((i/2)a ^ 



9r~^ /■+00 9^-1 „d/2 

^ e-«'^/"-'d. = P^-^r(d/a-l) 



(2^)'^ r((i/2)a io (2^)^^ r((i/2)a 

Appendix. 

Proof, (of Fact 12. 2p Let us first transform (jl3p using ([7]) 

/(t) := i-^/^-i (p,, Ao) = t-^/"-! / pj(x)Ao(dx) = / pi(t-i/"x)Ao(dx) 



Let us denote F{t) = Ao(i?(t)). The density pt of the symmetric a-stable process is rotationally 
symmetric hence we may denote p(|x|) := p{x). It is also know that it is non-increasing. We also 
denote Oj^, = ^t^^"- It is straightforward to check that 



n 

oo 



(46) t"' J^Pi {{i + l)/n) {F {at,i+i) - F {at,i)) < f{t) < t-^J^P^ (^/^) " ^ 



i=0 i=0 



Denote the right most part of the last expression by Dt^n- It is easy to check that t ^ {F (at,i+i) — F {at-, 
n{K) ((^)" - (^)") and (F(at,i+i) - F {at^i)) \ < Using the dominated Lebesgue con- ' 

vergence theorem we have 



i=0 

We use inequahties f (^)°~^ < (^)° - (^)" < f (^)°"^ for a > 1 (for < a < 1 one have to 
reverse inequahty signs) in order to prove 

r+oo 

hm Dn = n[A)a pi{x)x°'-^ dx . 

Using the same technique we proceed with "the left inequality" in (|46p . □ 

Acknowledgements. The author would like to thank Prof. T. Bojdecki for introducing him to the 
problem and Dr. A. Talarczyk for fruitful discussions. 
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